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We apply the non-local chiral quark model to study vector and axial 
pion-to-photon transition amplitudes that are needed as a nonperturbative 
input to estimate the cross section of pion annihilation into the real and 
virtual photon. We use a simple form of the non-locality that allows to per- 
form all calculations in the Minkowski space and guaranties polynomiality 
of the TDA's. We note only residual dependence on the precise form of 
the cut-off function, however vector TDA that is symmetric in skewedness 
parameter in the local quark model is no longer symmetric in the non-local 
case. We calculate also the transition form-factors and compare them with 
existing experimental parametrizations. 

1. Introduction 

Exclusive processes involving hadrons factorize in the Bjorken limit into 
a hard cross section and a soft hadronic matrix element that cannot be calcu- 
lated in perturbative QCD. Those matrix elements encode nonperturbative 
information on the hadronic structure. Recently Pire and Szymanowski 
[l] introduced new objects of this type that describe pion-to-photon (7127) 
transition in the presence of the qq operator that in the following will be 
denoted by T. Depending on the nature of T one can define vector or ax- 
ial transition distribution amplitudes (VTDA or ATDA respectively). In 
their original work Pire an Szymanowski discussed hadron-antihadron scat- 
tering process HH 7*7 where the virtual photon supplies the hard scale 
allowing for perturbative treatment, whereas the other photon is on mass- 
shell. As the simplest case, to avoid complications with spin or multiquark 
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bound states, one may take pions as initial hadrons. Although experimen- 
tally difficult to access [2], vrvr scattering is of particular theoretical interest, 
since pions are Goldstone bosons of broken SU(2) chiral symmetry and their 
properties are to large extent determined by the symmetry (breaking) alone 
rather than by the complex phenomenon of confinement. 

Indeed, there exists in the literature a variety of chiral models which 
involve both constituent quarks and pion degrees of freedom. In Ref.[3] 
Tiburzi used simple constituent quark model to discuss properties of the 7r27 
TDA's. Similar model with Pauli-Villars regularization has been recently 
used by Courtoy and Noguera [1] to calculate 7r27 TDA's for different ranges 
of kinematical variables. Finally Ruiz Arriola and Broniowski [5] employed 
the Spectral Quark Model (SQM) for the same purpose. 

One of the important ingredients of the low energy models is regulariza- 
tion. Even though 7r27 TDA's are formally finite, regularization cannot be 
simply dropped out, since it defines the scale of applicability above which 
the models do not apply. In this paper we calculate 7r27 TDA's in the 
semibosonized Nambu-Jona-Lasinio model known also as the Chiral Quark 
Model (xQM) with non-local regulator. This model has been previously 
used to calculate pion [6], pion and kaon [7] distribution amplitudes (DA) 
and generalized parton distributions (GPD) together with two-pion distri- 
bution amplitudes [8] . Direct comparison of local and non-local versions of 
the model allows to determine the influence of the non-local regulator. In 
most cases rather sharp curves obtained within the local model are smoothed 
down; also the endpoint behavior of various distributions is made contin- 
uous. This phenomenon is at best illustrated by the example of the pion 
distribution amplitude which in the local model is constant over the whole 
support, whereas the non-local regulator forces it to vanish in the endpoints 
[9l [To]. It is therefore of interest to investigate the role of the non-local 
regulator for the 7r27 TDA's introduced above. 

One has to remember that xQM, although devised to describe chiral 
physics of Goldstone bosons, has been widely used to incorporate baryons 
as chiral solitons both in local (for review see e.g. Ref. and non-local 
[12| cases. Generally the results of these studies show that the soliton ceases 
to exist for too small constituent quark mass M. The critical value of M 
depends on the details of the given model, however it is of the order of 300 
MeV or a bit less. Typical values of M that fit well the hyperon spectrum 
may be as high as 420 MeV [13]. In the present paper we adopt two distinct 
values of M: 350 MeV and 225 MeV. The latter, as we shah se in SectEl 
fits well the slope of the 7r27 transition form-factor. It is, however, excluded 
if one wants to describe baryons as chiral solitons. 

Our results can be in short summarized as follows: it seems that 7r27 
TDA's are quite robust as far as different regularization schemes (including 



tda 'app 'vl printed on January 30, 2009 



3 



no cut-off at all) are concerned. On the one hand this is a welcome feature 
for phenomenology, one can use them with large degree of confidence. On 
the other hand they cannot be used to distinguish between different models. 
However, as we shall discuss in Sect. 15.21 the ^-symmetry of the VTDA is 
no longer present in the non-local model, and the results for ,^ < differ 
more from the results of the local model, than the ones for ^ > (where 
^ is the skewedness parameter to be defined in Sect. U]). We shall also see 
that the normalizations of VTDA and ATDA which are equal in the local 
model become different in the non-local case. 

In the next section we shall discuss different chiral quark models existing 
in the literature in the context of 7r27 TDA's. Subsequently in Sect. [3] we 
give a short overview of the non-local model used in the present paper. We 
shall work in a symmetric kinematics that is introduced in Sect. H] together 
with the definitions of the TDA's in question. Calculations and results are 
presented in Sect. O We summarize and give our conclusions in Sect. [6l 
Technical details can be found in Appendices [XHUl 

2. Quark models and the transition amplitudes 

In order to estimate transition amplitudes one may try to construct 
phenomenological ansatze that satisfy general conditions such as gauge in- 
variance and anomaly structure, Lorentz invariance that in our case is equiv- 
alent to polynomiality, etc.. For special limiting cases these ansatze reduce 
to known form factors or structure functions (see e.g. Ref . [14] . or in the 
context of tt2j TDA's Ref. [3]). Alternatively one may try to resort to some 
kind of nonperturbative calculations. Usually a good staring point is a 
constituent quark model where all nonperturbative effects are parametrized 
in terms of a constituent quark mass and a chirally invariant meson-quark 
coupling (Local Chiral Quark Model - L^QM). As quite useful first ap- 
proximation such model has been used in Refs.[^ [TU] to describe pion light 
cone distribution amplitude or two pion distribution amplitudes [15]. The 
results of the LxQM seem rather trivial. For example pion light cone dis- 
tribution amplitude is constant and does not vanish in the endpoints [9l[6], 
similarly the isoscalar skewed pion distribution amplitude is just a super- 
position of the Q functions [6]. While LxQM satisfies Ward identities, it 
violates Lorentz invariance by the necessary transverse UV cut-off. 

Let us stress that the UV cut-off is not merely a regulator, on the con- 
trary, it defines the applicability domain of the model that is clearly devised 
as low energy approximation to QCD. It should be therefore applied also 
to the quantities that are formally finite in this limit. However, as we shall 
shortly explain in some more detail, the explicit UV cut-off results in a vi- 
olation of the polynomiality which is essential for the transition amplitudes 
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that are discussed in this paper. Clearly a more sophisticated regulator is 
needed. 

One might expect that a more sophisticated regularization which pre- 
serves Lorentz invariance would result in a more realistic shape of the pion 
distribution amplitude. This is, however, not necessarily the case. Gauge 
invariant regulator is provided for example by the Spectral Quark Model 
(SQM) Refs.[16l [5], where the constituent quark mass M is traded for a 
spectral parameter iv, and all physical quantities are given in terms of inte- 
grals over du> with some a priori unknown spectral density p{uj). Spectral 
density p{io) must satisfy a number of relations that follow from the QCD 
Ward identities; explicit realizations of the model with explicit form of p{uj) 
are also known [IB]. Although at first sight theoretically attractive, SQM 
yields phenomenological results that are very similar to the naive LxQM 
described above. 

Similarly in Ref.[T7] a Pauli-Villars regularized Nambu-Jona-Lasinio 
model was applied to calculate both pion distribution amplitude (DA) and 
parton distributions in the pseudoscalar mesons. Again pion DA at the in- 
put scale is given as a step function and does not vanish in the endpoints. 
Evidently a form factor in the quark-pion vertex is needed to "soften" the 
shape of the (generalized) distribution amplitudes. 

An attractive and simple way out is provided by the Non-Local Chiral 
Quark Model (NLxQM) where the quark-pion coupling is given in terms of a 
momentum dependent constituent quark mass M (p). For small p, M{p) 
M ~ 350 MeV, whereas for p oo, M{p) 0. This suppression of large 
momenta (remember that the constituent quark mass M{p) acts not only 
as a mass parameter in the propagators, but also - more importantly - as a 
quark-meson coupling) is enough to make the pion distribution amplitude 
vanish in the endpoints [6]. 

Momentum dependent constituent mass M[p) preserves polynomiality, 
however it violates QCD Ward identities. The latter can be easily under- 
stood by the following example. Consider Dirac equation with a momentum 
dependent mass M{p): 

[i> - M{p)) u{p) = (1) 
and the electromagnetic current 

f =u{p + qh^u{p). (2) 

Naive current conservation reads 



= u{p + q) [{^ + ^) - u{p) 

= [M{p + q)- M (p)] u{p + q)u{p) / 0. 



(3) 
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Modifications of tlie electromagnetic current tliat make tlie current 
conserved liave been proposed in Refs.[T8] -[27j- The discussion of the low 
energy theorems in the context of the instanton model leading to the mo- 
mentum dependent constituent quark mass with special emphasis on axial 
anomaly can be found in Ref . |28j . Although these modifications, supplied by 
an appropriate requirements of the absence of kinematical singularities [19] 
fix the longitudinal part of the pertinent non-local vertices, the transverse 
part remains undetermined. When one wishes to consider electromagnetic 
processes, it is necessary to assume some model for the transverse part. 
Therefore for the purpose of the present work we do not consider such mod- 
ifications, although we acknowledge the fact that such a study is certainly 
required. Below we give an argument in favor of such a procedure that 
follows from the parametrical dependence of the current nonconservation in 
the instanton model of the QCD vacuum. 

A non-local model leading to ([T]) can be "derived" from QCD in the 
instanton model of the QCD vacuum In this model the vacuum is 

filled with interacting instantons that stabilize in a configuration where 
the mean instanton radius p ~ 1/(600 MeV), whereas the typical instaton 
separation R ~ 1/(300 MeV). The instanton packing fraction [p/R)^ is a 
dynamical small parameter of the model. The model allows to calculate 
M{p) in Euclidean space. Setting M{p) = MF^{p) we have 



where z = pp/2 and M is the constituent mass at zero momentum. There- 
fore (schematically) 



Equation ([5]) maybe viewed as an expansion in the inverse momentum 
Qinst. = 2/;0 corresponding to the typical instanton size. Hence for small 
momentum transfers (and this is certainly the domain of the present model) 
the nonconservation of the vector current is parametrically small in the in- 
verse instanton size Qinst.- Therefore in the following we shall use local 
currents, such as ([2]), rather than the non-local extensions, allowing for the 
violation of Ward identities at the level of q/Qinst.- The price we pay for 
that is wrong normalization of the pertinent form factors, since it is fixed by 
the axial anomaly. However, the dependence on the kinematical variables 
is almost identical as in the models that preserve Ward identities. We shall 
come back to this point in Sect {61 Finally let us stress that despite the fact 
that our currents do not satisfy Ward identities the amplitudes we calculate 



Fi„st.(p) = 22 [Io(z)Ki{z) - h{z)Ko{z)] - 2h{z)Ki{z) 



(4) 



M{p + q)- M{p) = M [F\{p + q)p) - F\pp)\ 



2 dz 



(5) 
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are gauge invariant, in the sense that they vanish when contracted with 
on-shell photon momentum. 

3. Non-Local Chiral Quark Model 

In order to provide non-local nonperturbative regulator we employ semi- 
bosonized Nambu-Jona-Lasinio model defined by the following action de- 
scribing quark interaction with an external meson field U \29\ [30] : 

/rl^kH^l - 
-^^{k)F{k)U^^ {k - l)F{l)m (6) 

and U^''{x) can be expanded in terms of the pion fields: 

W^^{x) = 1 + ±-^\^^\x) - ^vr^(x)7r^(x) + . . . (7) 

M is a constituent quark mass of the order of 350 MeV and F{h) is a 
momentum dependent function such that -F(O) = 1 and Fik"^ oo) 0. 
In what follows, for comparison, we will consider also M = 225 MeV 

Note that ([6]) provides both momentum dependent mass of the quark 
fields and the non-local quark-meson coupling. Pions act at this stage only 
as auxiliary fields being - by equations of motion - objects composed from 
quark-antiquark fields. Kinetic term for pions appears only after integrating 
out the quark fields [301 EI] and the proper normalization is obtained by an 
appropriate choice of the cut-off function F{k). Here, following Refs.[9l[6] 
we wish to perform all calculations in the Minkowski space. To this end we 
choose: 

^(^-) = {w^^ («) 

which reproduces reasonably well @ for A;^ < 0. Numerical values of A„ 
for different choices of M are given in Table I of Ref . [6] and Table [1] of the 
present paper. Scale A„ should not be confused with the typical momentum 
scale Qq (which for the original shape of F{k) given by Eq.(|3|) is equal to 
Qinst = 2//o), that can be defined as the value of the momentum for which 
F{Qo) = const., say 1/2. Then 

Qo(n) = (9) 

and does not exceed 2 GeV for the highest values of A^. 

Ansatz ([8|), apart from being close to the instanton motivated function 
([5]), is very practical for calculations in the Minkowski space. Indeed, it 
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introduces a number of complex poles in the complex momentum plane, 
that can be analytically integrated over in the light cone coordinates. Light 
cone coordinates are defined by two null vectors: h = (1,0,0, 1) and n = 
(1, 0, 0, —1). In this kinematical frame any four vector v can be decomposed 
as: _^ 

v^' = ^n^^ + ^n^^ + (10) 

with u"*" = n • f , v~ = n ■ V and the scalar product of two four vectors 
reads: 

1 , _ 1 _ , ^ ^ 

V ■ W = -V W +2"^ W —VT-WT- (11) 

Therefore 

k'^ -Al = k-k+ -k^-Al. (12) 

The integration measure in the light-cone coordinates takes the following 
form: 

d'^k = dk-^ dk~d^kT/2. (13) 

Looking at (jl2p we see that ([8]) generates a n— th degree pole in the k~ 
plane that can be easily integrated over. The details can be found in Ref.[6] 
and in SectiH 



4. Definitions and kinematics 



We use the definitions of pion TDA's from [3] (our definitions include 
additional i phase factor) 



2^' 



X {l{P2,e) 
ie 



A 



d I n I "f^u I —n 



A 



vr+ (Pi) 



(14) 



dX , 



27r 



X (7(^2,e) 
e 



-n 



P^{q-e*) A{X,C,t) + ... 



'2^/2F.p+-^^ ^.^^.^,s,^. , ... (15) 

where V {X, ^, t) and A {X, ^, t) are vector and axial TDA respectively (we 
use i^TT = 93 MeV). We shall use the following symmetric parametrization 
of momenta: 



Pi 



:i-e)Y^M + (i+e) 



T 



2p+ ' 



(16) 
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where px = 0, with 

p, = ^ (Pi + P2)^ = ^h, + ^n, (17) 
Here ^ denotes skewedness. Momentum transfer reads: 

9m = (^^2 - Pi)^ = -ep+n^ + e^^M + qjl- (18) 

Note that is related to p'^. Indeed, using the on mass-shell condition (for 
= 0): 

we arrive at: 

Considering the momentum transfer squared: 

= t = -4fp^ - ql or - 4^ V = (t + , (20) 

and inverting (j20p . we get an important constraint: 

= -(1 - ^2)t > ^ -1<C<1. (21) 

One has to note, that the lower limit for ^ depends on the order of limits at 
t = and rriTr / [4]. We have avoided this ambiguity by keeping nin = 
from the very beginning. 

Photon polarization vector satisfies e* ■ P2 = 0. Dots in Eq. (jl5p stand 
for parts which are structure independent and thus irrelevant in our consid- 
erations [Sin]. In general TDA's do not posses any symmetry properties in 
^ and X - in contrary to GPD's [3]. 

One can define flavor diagonal VTDA's I4 i^-, t) and {X, ^, t) by re- 
placing in definition vr"'" by 7r° and taking operators u {—^nj ^^u 
and d ( — ^n) 7'^(i (^n) respectively. On the other hand matrix element 
{'j{P2,e) (0) (0)| 7'"'' (-Pi)), where -0 are now iso-doublets, is param- 
eterized by the pion-photon transition form factor controlling 7*7 — > 
processes. Therefore one can derive the sum rule pQ relating this form fac- 
tor to Vu {X, ^, t) and {X, ^, t). Its normalization is fixed by axial anomaly 
connected with the Ward identity relating matrix elements for transition of 
axial current to two photons with a similar matrix elements of a pseudo- 
scalar current. This normalization together with the conventions of Eq. (jl4p 
gives the normalization condition for the VTDA: 

1 

ldXViX,^,t = 0) = ^{Qu + Qd) = ^, (22) 

-1 
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which is independent of M and of ^. The latter is related to the polynomial- 
ity which states that the n—ih moment of the TDA's in X is a polynomial in 
^ of degree not higher than n. Normalization (|22p is automatically satisfied 
in the local chiral quark model which in the chiral limit {i.e. for jn^r = 0) 
gives the normalization of the ATDA of Eq. (jl5p equal to the one of the 
VTDA: 

1 

jdXA{X,i,t = Q) = ^{Qu-Qd) = i^. (23) 



Note, however, that (j23p is not fixed by the anomaly. This will be important 
in the non-local model where the two normalizations are not equal any more. 

Moreover, we have the following sum rules relating vector and axial- 
vector form factors with the relevant TDA's 

fdXD (X, e, t) = ^^^Fd (t) = 2V2F^F^ (t) , (24) 

-1 

where D stands for V or A. In the chiral limit, which is considered in this 
paper, we define Fy (t) = Fy (t) /m^r and similarly for F4. 

5. Transition distribution amplitudes in the chiral quark model 

Using effective action ([6]) we obtain the following expressions for the 
matrix elements (1141) and (1151): 



d ( --n] T^'u i-n 



vr+ (Pi) 



= {Q.MI" + Q.Mf ) e:, (25) 

where F'^ is either 7^ or 7^75, is number of colors, Qu and are charges 
of quarks u and d respectively. Two amplitudes Aii and A^2 depicted in 
Fig(T]are defined as (we omit ie prescription in fermion propagators): 

= ( -^,5{k+-{X-l)p+)F{k)F{k + P,)^ (26) 
J (27r) 

Tr i i r'' ^ 7^ i 7^ 

^2-M{k + P2) jk+ /Pi-M{k + Pi) /c-M{k)' 

= f 4^Hk^-iX + l)p+)F{k)F{k-Pi)x (27) 
J (27r) 

Tr I F'' 7^" 7^ 1 . 

1 /:-/'! - M {k - Pi) ^- /P2 - M (/c - P2) M (k) J 
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r 



k=k+P, 



u 



+ 




k 

r 



k 



Fig. 1. Feynman diagrams for J\4i and A42. Traces should be taken opposite to 
the momentum flow denoted by arrows. Time flows from right to left. 



First we calculate TDA's neglecting the mass dependence upon p, i.e. 
we set M (p) = M. It was already done within a very similar model in 
Ref. [1], where the calculations were performed in the Minkowski space 
with Pauli-Villars regularization. In another approach partially discussed 
in Refs.[9]-[8] one uses the light cone coordinates (jlOh with the integration 
measure given by Eq. (jl3p . performing first integration over dk~ and then 
over dk^ with the transverse cut-off A'^. The latter is chosen to normalize 
appropriately F.,^ or alternatively to normalize the pion distribution am- 
plitude to 1. Introducing transverse cut-off in the integrals defining TDA's 
would violate Lorentz invariance and, as a consequence, polynomiality. This 
can be nicely illustrated by considering the integral 



which is related to the zeroth moment of the vector TDA and expanding it 
for small t: 



5.1. Local case 




(28) 
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8(27r)2 y (4 + M2)2 '47 ^ (4 + M2) 

Now, in order to make dependent term vanish 

2kl-M'^ _ A2 

"^^'^(4 + M2)4 - "(M2+A2)3 (30) 


we have to choose A2 = oo, what leads to: 

^= 8(27r)2M2|^ + T2M2 +■■■} ^^^^ 

where . . . denote higher powers of t. 

So in order to preserve polynomiahty we have to work with an infinite 
transverse cut-off. In this case it is very useful to switch to Euclidean space 
and use Schwinger representation for scalar propagators, following Ref.[5]. 
This allows to obtain analytical results in a very simple way. We shall be 
using TDA's calculated in the local model as a reference when discussing 
the results in the non-local case. 

Calculating traces and combining definitions , (fT5]) with ([25|) we get 



V (X, e, t) = mM^N,p+ (Qd/Ci + QuK,2) , (32) 

A (X, e, t) = '—iQM^N,p+ {QdJi + QuJ2) , (33) 

q ■ e* 



where 



(2vr) ({k±Pif-MA{{k±P2f-MA{k-^-M^) 



(2vr) ((A;±Pi)2-M2) ((A:±P2)^-M2) (fc2-M2) 



with upper signs referring to subscript "1" and lower signs to "2". 
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Taking the same steps as in [5], we get 



{47rfp+Io '^^lo 



X 5(y(l+0p+ + 2(l-0P+±(^Tl)) 



1 



(36) 



M2 - yzt 



J\,2 




X 5(y(l+0p+ + 2(l-0P+±(^Tl)) 



\-2y 



(37) 



M2 - yzt ■ 



While obtaining the second equation we used n • e* = (in the Ught-cone 
gauge), what made that expression finite. Let us notice that from these 
formulae it is obvious that TDA's satisfy polynomiality condition. Simple 
integration over Feynman parameters leads to the final result, which we 
quote in Appendix [XI 

The transition form factor we get in the local model recovers the nor- 
malization required by the axial anomaly. Its analytical form - calculated 
long time ago in Ref. [32j in more general kinematics - is given in Appendix 



In this section we take full momentum dependence of the constituent 
quark mass in integrals (j27p and (j28p . The calculations will be done in 
the Minkowski space. We use method of evaluating the contour integrals 
developed in Refs.[Bl[S]. 

We start from VTDA. When evaluating the trace we find not only struc- 
tures proportional to s^^'^^ e\,'Paq(i^ which are needed, but e^^°^^ e\,naP\(i and 
e^'^°'^ £*naP2i3 as well. However they vanish when contracted with n^, thus 
they are unphysical in the sense that they do not give contribution to the 
observables. We obtain expressions of the following form 



El 



5.2. Non-local case 




where 



D (p) =p'^ -M"^ (p) + ie. 



(39) 
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Fuctions Ai^2, -^1,2) (^1,2 depend on X, ^, t and integration variable kx- 
Their explicit form is given in Appendix [Cl 

Next, we have to take mass dependence on momentum given by ([8]). We 
choose Aji for given n in such a way that pion DA calculated in the present 
model is normalized to unity. Pertinent values of A„ are listed in Table [H 



n 


12 3 5 


M = 350 MeV 
M = 225 MeV 


1156 1727 2155 2819 
2121 3125 3880 5060 



Tabic 1. Values (in MeV) of A in function of n for M = 350 and 225 McV. 

Introducing dimensionless variables k = k/A, Pi = -Pi /A, P2 = -P2/A, 
r = M /A and using 

= (k±A,2)^- 1 + ie, (40) 

n3 = K^-l + ie, (41) 



we get 



1 r d?KTdK-dK+ 6 (k+ - (X =f 1) p+) 
2MA3 J J^"^ G(n±)G(n±)G(n3) 

+B,,2 {ufriutriuir 

+Ci,2(nf)^"(4)'"(4)"}, (42) 



where G (u) = + u^"- — r^. Polynomial G (u) can be alternatively 

written in a factorized form G (u) = H^i^^ {u — Zi) , where Zi are roots of 
equation G («) = and can be found numerically. Note, that if r = 
(i.e. A — > 00) we have 4n degenerate solutions equal to zero and one equal 
to —1. If A becomes finite the degeneracy is lifted and we have 4n + 1 
solutions which in general are complex. Integration over dK~ has to be 
done by the residue theorem, thus we have to find the poles in complex 
plane. However, because of the imaginary part of z/s, the poles can cross 
the standard integration contour. This may result in non- vanishing of the 
TDA's in unphysical regions. To avoid this, the integration contour has 
to be modified. Detailed discussion of these problems is given in [6] and 



14 



TDA'APP'VI PRINTED ON JANUARY 30, 2009 



Appendix After performing the contour integrals we get: 

4n+l „ ,2 



i,j,k=l 



» Pi, 2 ri2 s-12 , D "'i.a ^1,2 ci,2 , ^ di 2 ei,2 91,2 
iai,2Zi - Zj + ^1,2) (7l,2^;i - Zk + Pl,2) 

where 

4n+l 

/.= n (^^-^^•)"' (44) 





Fig. 2. Vector TDA's in the non-local model: a) for n = 1 and ^ = 0.2 (dashed), 
5 = 0.5 (solid), ^ = 0.8 (dash-dot); b) for n = 1 and S, = -0.2 (dashed), ^ = -0.5 
(solid), ^ = —0.8 (dash-dot). Comparison of the local model: c) for ^ = 0.5 
(dashed) with the non-local model for n = 1 (solid) and n = 5 (dash-dot); d) for 
5 = —0.5 (dashed) with the non-local model for n = 1 (solid) and n ~ 5 (dash-dot). 
All plots are made for constituent quark mass M = 350 MeV and t ~ —0.1 GeV^. 
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Fig. 3. Vector TDA's in the non-local model: a) for n = 1 and ^ = 0.2 (dashed), 
f = 0.5 (solid), f = 0.8 (dash-dot); b) forn = 1 and ^ = -0.2 (dashed), ^ = -0.5 
(solid), = —0.8 (dash-dot). Comparison of the local model: c) for = 0.5 
(dashed) with the non-local model for n = 1 (sohd) and n = 5 (dash-dot); d) for 
f — —0.5 (dashed) with the non-local model for n ~ 1 (solid) and n = 5 (dash-dot). 
All plots are made for constituent quark mass M — 225 MeV and t = —0.1 GcV^. 



Here powers of Zj-s (denoted by Latin characters), as well as the explicit 
form of the functions ai^2, Pi,2, 7i,2, Pi, 2 depend on the region of X (see 
Appendix E]). It should be pointed out that expressions denoted by Greek 
characters contain second power of kt, while functions ^i,2i ^i,2) ^^1,2 are 
of first order in kt- Therefore the integration over (Pkt = nTdnxdOT is 
finite. Integral over (IOt can be done analytically by integration over a unit 
circle, while integral dKT can be performed numerically. 

In the non-local model we cannot recover the required normalization 
(j22p for finite A. This can be understood as a consequence of the regular- 
ization that does not respect axial anomaly. Therefore we impose the proper 
normalization (j22p by multiplying the VTDA by a suitable correction factor 
Ny as given in Table [2l 
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Results for VXD A are shown in Figs 12] and [3] for constituent masses 
M = 350 and 225 MeV respectively. We observe that the "non-local" 
curves are less sharp than in the local case and that their shape depends 
slightly on n. Also the maxima are shifted from X = it^, where they were 
placed in local case. The non-local model has the feature that VTDA's are 
no more ^-symmetric and for ^ < it gives results that are more peaked 
than in the local case and with the middle zero in a different position. The 
deviation from the local model is stronger for larger constituent masses. 





Fig. 4. Axial TDA's in the non-local model: a) for n = 1 and ^ = 0.2 (dashed), 
C = 0.5 (solid), ^ = 0.8 (dash-dot); b) for n = 1 and ^ = -0.2 (dashed), ^ = -0.5 
(solid), ^ = —0.8 (dash-dot). Comparison of the local model: c) for ^ — 0.5 
(dashed) with the non-local model for n = 1 (solid) and 7i = 5 (dash-dot); d) for 
5 = —0.5 (dashed) with the non-local model for n = 1 (solid) and n = 5 (dash-dot). 
All plots are made for constituent quark mass AI ~ 350 MeV and t = —0.1 GeV^. 

In the case of axial TDA's the algebraical steps are the same. The 
complication is that after evaluating the trace we have to retain only terms 
proportional to Pg* ' since all other terms are structure independent 
or gauge artifacts. General expression is similar to (j44p . but now functions 



tda 'app 'vl printed on January 30, 2009 



17 



A, B, C contain also second power of kt- However, the integration over (Ikt 
is finite because of the property: 



4n+l 



i=l 

(see again [6] and Appendix ICl). 



if m < 4n 

1 if m = An 



(45) 
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Fig. 5. Axial TDA's in the non-local model: a) for n = 1 and ^ = 0.2 (dashed), 
C = 0.5 (solid), £, = 0.8 (dash-dot); b) for n = 1 and ^ = -0.2 (dash-dot), ^ = -0.5 
(solid), ^ = —0.8 (dashed). Comparison of the local model: c) for ^ = 0.5 (dashed) 
with the non-local model for n = 1 (solid) and ?? = 5 (dash-dot); d) for ^ = —0.5 
(dashed) with the non-local model for n = 1 (solid) and n = 5 (dash-dot). All 
plots are made for constituent quark mass M = 225 MeV and t = —0.1 GeV^. 

Numerical results for the ATDA's are shown in FigsU] and [5] for con- 
stituent masses M = 350 and 225 MeV respectively. W see again that to a 
good accuracy all models (local and fully non-local one for different n) give 
the same results both for positive and for negative ^ (note, however, small 
shift of the minima in the non-local case) . All curves were normalized as in 
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the local case according to Eq. (l23p . multiplying the calculated distribution 
by correction factors listed in Table [3l Again the deviation from the local 
model is stronger for larger constituent masses. 

We have checked numerically that TDA's in the non-local model satisfy 
polynomiality condition for first three moments. However, contrary to the 
local case, even moments of VTDA are nonzero. 

6. Summary and discussion 

In the present paper we have employed chiral quark model with momen- 
tum dependent constituent quark mass to calculate pion-to-photon transi- 
tion distribution amplitudes. Before we briefly summarize our results let 
us discuss the main features of the model. We have chosen momentum 
dependence in the simple form given in Eq.dS]) which for Euclidean mo- 
menta resembles M(p) obtained within the instanton model of the QCD 
vacuum. This form of M{p) allows to perform all integrations directly in 
the Minkowski space and has been previously applied to calculate pion [6] 
and kaon [25] distribution amplitudes and two pion distributions as well as 
generalized parton distribution of the pion [8] . The main technicality that 
we wish to mention, consists in the proper choice of the integration contour 
in the loop momentum k~ which is discussed in Appendix O and can be 
also found in Ref.[6]. 

Proper choice of the integration contour guaranties that the TDA's are 
real and have proper support in kinematical variables X and ^ defined in 
Sect m and satisfy polynomiality. 

Throughout this paper we have used momentum dependent mass that 
acts as an UV cutoff and, at the same time, as the quark form factor within 
the pion. The latter is very important for making the pion DA vanish in 
the endpoints [6l |9]. Indeed, similar calculation of the photon DA yields 
distribution that is discontinuous in the endpoints [lOj reflecting the point- 
like nature of the quark-photon coupling. 

Although we have used momentum dependent mass, we have not modi- 
fied currents accordingly [l8]-[27j, and as a consequence our model violates 
QCD Ward identities. This violation is however "mild" as it occurs at the 
level qp/2 where p is the mean instanton radius (0). Nevertheless viola- 
tion of the axial Ward identity results in the wrong normalization of the 
vector TDA which is fixed by the axial anomaly. In order to get over this 
deficiency we have simply corrected normalization of VTDA to the value 
obtained in the local model which does satisfy axial Ward identity. The 
correction factors (model results should be multiplied by Ny or to ob- 
tain normalization of Eqs. (j22l23p ) are given in Tables [2] and [3] and may seem 
large. Local limit can be obtained by pushing artificially A — > oo in (l8|). 
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Obviously the same procedure applied to the pion DA would yield pion DA 
constant, but with an infinite norm. In that case correction factor would be 
infinite (modulo some regularization such as a transverse cutoff for exam- 
ple). On this scale correction factors of the order of 1.5 are not excessively 
large. 



M 


n = 1 n = 2 n = 3 n = 5 


225 MeV 
350 MeV 


1.151 1.148 1.147 1.146 
1.487 1.490 1.490 1.491 



Table 2. Correction factors Ny for VTDA for AI = 225 and 350 MeV and different 
values of n. 



M 


n = 1 n = 2 n = 3 n = 5 


225 MeV 
350 MeV 


1.083 1.081 1.081 1.080 
1.217 1.219 1.219 1.219 



Table 3. Correction factors Na for ATDA for M = 225 and 350 MeV and different 
values of n 

Despite the fact that Ward identities are not satisfied our amplitudes are 
gauge invariant, i.e. they vanish when contracted with photon momentum. 

Normalization of axial TDA's is not fixed by the anomaly. However, in 
order to compare them with local model we used correction factors fixing 
normalization given by (j23p . This normalization overshoots experiment by 
approximately a factor of 2. Such a large mismatch is common to local 
quark models [5]. The correction procedure used in this paper to maintain 
normalizations (j22p and (j23p is to large extent arbitrary. Taking normal- 
izations as they come out (i.e. without correction factors Nv^a) would shift 
the axial transition form factor at t = towards the experimental value. 
At the same time VTDA would loose correct normalization. The latter, 
however, should be attributed to the violation of the Ward identities in 
the present version of the model and it is the violation of the axial anomaly 
which is responsible for the wrong normalization in the vector case. Clearly, 
only a complete calculation with the non-local currents might resolve this 
discrepancy. 

Transition form factors are defined in Eq. (j24p . We show them in Figl6j 
All calculations were performed for the constituent quark mass M = 350 
and M = 225 MeV. We find that in the case of the non-local model the 
transition form factor is more dumped than the one calculated in the local 
version. However, for realistic constituent quark mass M = 350 MeV it 
still falls off much slower than the experimental curve parameterized by the 
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function [33] 

where Mq = 776 MeV. We could get good description of experimental data 
for much lower values of constituent quark mass parameter M. This can be 
easily understood from the approximate formula (j3ip where the slope reads 

= 12M^ M = 225 MeV. 

For that reason we have used M = 225 MeV for our calculations although 
the reasonable values of the constituent masses that have been used in in 
the literature - as explained in Sect. [T]- lie above 300 MeV. 



0.28 -, 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 
- 1 [GeV'] 



Fig. 6. Comparison of the transition form factors obtained in various versions of 
the chiral quark model. Solid line represents experimental fit of Eg. ((46)) . Dashed 
lines correspond to local model, dash-dotted lines to the non-local one with n — 1 
(there is almost no n dependence). Two upper curves correspond to M ~ 350 
MeV, two lower ones to M = 225 MeV. 

Our calculations were performed in the symmetric kinematics defined 
in Sect iH and can be directly compared with Ref.[l]. There is qualitative 
agreement between the local version of the present model and the one of 
Ref.[l] (up to an overall sign of ^ for the axial case, see Ref . [34j ). In order 
to make comparison with [5] we have repeated their calculations in our 
kinematics. The results are essentially identical to the local version of the 
present model, provided we take small constituent mass. This is illustrated 
in Figs. [7] where we plot vector and axial TDA's for ^ = ±0.5 and t = —0.3 
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GeV^ in SQM and local version of the present model for M = 225 and 350 
MeV. 
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Fig. 7. Comparison of vector (a) and axial (b,c) TDA's obtained in the SQM of 
Ref. [5] (solid) with the local version of the present model for M — 225 MeV 
(dashed), M = 350 MeV (dash-dot) and £_ = ±0.5 (a), 0.5 (b) and -0.5 (c). Plots 
are made for t = —0.3 MeV. 

We see that the pion-to-photon transition amplitudes defined in Eqs. (fH|) 
and ([15]) are quite robust. Basically their shapes do not depend on the 
specific model and on the regularization used. Nevertheless some small 
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differences between local and non-local models can be observed. The most 
prominent is the violation of the ^—symmetry present in the local case 
for the vector TDA which can be seen in Figsl2] and [3l For negative ^ 
the non-local model gives results that are more peaked than in the local 
case and with the middle zero in a different position. On the other hand 
ATDA's are very close to the local case. One may be therefore confident 
that the shape of the axial TDA is without doubts as shown in FigsjH and 
[5l however normalization is not certain and should be perhaps adjusted to 
the experimental value of the axial transition form factor. 

Prom the point of view of QCD the quantities we calculate depend on 
a nonperturbative scale Qq which, however, must not be confused neither 
with the constituent mass M nor with an auxiliary parameter A. For A;^ < 
the Ansatz ([8]) should imitate M{k) obtained from the instantons. And for 
the latter, as explained in the Introduction, Qq ~ 2/p = 1200 MeV. It is 
therefore natural to assume that Qq is of the order of 1 GeV irrespectively 
of M and A. The precise definition of Qq is only possible within QCD 
and in all effective models one can use only qualitative order of magnitude 
arguments to estimate Qq. Discussion of this point can be found in Ref.[6]. 
Once the nonperturbative scale Qq is fixed, our results should be evolved 
to the hard scale characterizing given experimental setup by means of the 
evolution equations discussed recently at length in Ref. |35j . This will be a 
subject of a separate study. 
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Appendix A 

Transition distribution amplitudes in the local model 

According to notation introduced in Sect. 15.11 we obtained the fol- 
lowing expressions for VTDA (we quote results for positive ^, due to the 
^ — symmetry) 



(47r) p^ 



e{x + i)e{^-x) 



sign (ai,2 - Ci,2) 



In 



(ai,2 - ^1,2) (61,2 + Ci^2 



(ai,2 + ^1,2) (&1,2 - Ci^2) 



VAB 



1,2 



ai,2 — Cio 



0'l,2 + Ci^2 



(A.l) 



where ^ = (1 + ^) i, Ci,2 = y^WJjA. and 

l^X 



ai,2 



2(1 + 



(A.2) 



^1,2 



B 



1,2 



2e(i + e) 

(X^ 1)2^2 



4^ 



(1-0 M\ 



(A.3) 



(A.4) 



We remind that upper signs refer to subscript "1" . For ATDA we have 
for positive ^ 



Ji,2 (e > 0) = {e{x + i)e{i-x) + ^^^^ ("^'^ - 



X In 



(47r) p+ 
(ai,2 - ^1,2) (fei,2 + (^1,2) 



(ai,2 + (^1,2) (^1,2 - Ci^2) 



A 



-In 



2./AB 
Aaj 2 - ^1,2 



1,2 



(IT-)P(T- .)a + -Pi.2)^ign(ai.2-Ci,2.^^ 



Aftf 2 - Bi,2 

0-1,2 — C*! 9 



1,2 



ai.2 + 



(A.5) 
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and for negative ^ 



X In 



(47r) p+ 



(ai,2 + Ci,2) (fel,2 + Ci,2) 



Aaf 2 - Bi^2 



2 - Bi,2 
«1,2 — C12 



ai,2 + (^12 



(A.6) 



We introduced above 



X - 1 



D2 



X + 1 



1+e ' i+c 

All the remaining notation is the same as in the VTDA case. 



(A.7) 



Appendix B 

Transition form factor in the local model 



We obtain the following expression for the pion-photon transition form 
factor 



-F'-n-'y (t) 



M2 1 



27r2F^ t 



Li2 ( — ) +Li2 ( — 
a+ J \a- 



where 



-4 -v-^ 



(B.l) 



(B.2) 



and Li2 (x) is the dilogarithm function, defined as Li2 (x) = — Jq dt. 



Appendix C 

TDA 's in non-local model 



In FiglS] integration contour on the k~ complex plane is schematically 
shown. As explained in more detail in main text, contour is chosen in such 
a way that poles form each group cannot cross it. 
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Fig. 8. Integration contour on the k~ complex plane for r = M/A ^ 0. When 
r ^ the poles from each group merge into one pole and contour becomes standard 
semicircle. 

In the case of VTDA we have the following general formula for the 
integrals ([38]) 

i ^"^"^ f cPk 

i,j,k=l ^ ' 

A Pi, 2 ri,2 SI, 2 I D f^i-a ^1,2 ci.2 , di^2 ei,2 gi,2 

{ai,2Zi - Zj + /3i,2) {ll,2Zi - Zk + pi,2) 

First consider ^ > case. For Ai^2, Bi^2, (^1,2 we have the following 
expressions: 

Ai2 = ^-{X^l)^—^-^ (C.2) 

(C.3) 

The explicit form of ai.2, /3i,2i 7i,2i Pi, 2 and £1,2 depends on the region 
of the support under consideration. We introduce 

,1,2 _ (XTl)(l-6t T • + 1 ,^ 
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p+ {X ^ 1) 
In every region we have: 



-1<X <-i 



i<x<i 



ei = 0, 62 



X + l 



x + i x-i 



/?2 = p+ (X + i) {l 

P2 = p+(x-e) {ui-ui) 



1 

ei,2 = =F 



ei = ^, 62 = 

^-^ x + i 
A = p+ (X - i) {ui - ui) 



pi =p+[X + { 



Powers of Zi-s in the numerator can also be different in 
They read 
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-1 < X < and ^ < X < 1 



Pi,2 = 3n, ri,2 = 4n, si^2 = n (C.20) 
«i,2 = 3n, hi,2 = 2n, ci^2 = 3n (C.21) 
di,2 = n, ei,2 = 4n, 51,2 = 3n (C.22) 



Pi 


= n, 


ri = 3n, 


Sl = 


An 


ai 


= 3n, 


bi = 2n, 


Cl = 


■ 3n 


di 


= 3n, 


ei = n, 


91 = 


An 


P2 


= n, 


r2 = 4:n, 


S2 = 


3n 


a2 


= 3n, 


62 = 2n, 


C2 = 


: 3n 


d2 


= 3n, 


62 = 4n, 


92 = 


= n 



(C.23) 
(C.24) 
(C.25) 
(C.26) 
(C.27) 
(C.28) 

For negative ^ we have to change from Eqs. (|C.8IIC.l"9|) only Eqs. ()C.12I - 
IC.lSp . Appropriate expressions read 

ei. = T^ (C.29) 

«i,2 = ^, 71,2 = ^ (C.30) 

/3i,2=p+(X±e)(^^^'-ni'2) (C.31) 

pi,2=p+(XTl)(n^'2-ni'2). (C.32) 

Also Eqs. (|(L23II(:;.28I) have to be replaced by 

pi,2 = 4n, ri,2 = n, si,2 = 3n (C.33) 

ai^2 = 2n, 61^2 = 3n, ci^2 = 3n (C.34) 

di2 = An, 61,2 = 3n, 51,2 = (C.35) 

In the case of ATDA the general formula for J7i,2 is the same as for 
VTDA, but now we have different expressions for Ai^2^ ^i,2i Ci,2- They are 

Ai^2 = 9^^ (C.36) 
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Bi,2 = -f^ - 9^ (C.37) 

Ci,2 = -2g^T2v^ + f^±l, (C.38) 

where 

/- = ±l + ^-^f^ (C.39) 

9^ = --iX±l)-^P^ (C.40) 



2(i-o(e+i)'tn 

+ 2A (2^2 +^-l)t{X±l)KT-qT 

+ {e-l)t[{C + l){X±lft-4AC4) |. (C.41) 

Powers of Zi-s are the same as in the vector case. 
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